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ON STABILITY OF THE IN-HUMAN HOST AND IN-MOSQUITO 
DYNAMICS OF MALARIA PARASITE 


MOHAMED A. SELEMANI!?, LIVINGSTONE S. LUBOOBI!? AND YAW NKANSAH-GYEKYE! 


ABSTRACT. Stability analysis of dynamical system is basic requirement for its applica- 
tion in real-life settings. However, investigation for local stability is simpler than that 
for global stability, though the latter is more preferable. In this study, we perform sta- 
bility analysis of mathematical model for in-human host and in-mosquito dynamics of 
malaria parasites, and establish the existence of two types of equilibrium: malaria-free 
equilibrium (MFE) and malaria-infection equilibrium (MIE). Using linearization of sys- 
tem, MFE is proved to be locally asymptotically stable. By Metzler matrix theory, the 
MFE is reported to be globally asymptotically stable provided Ro < 1. By applying 


Lyapunov functional method and LaSalle’s invariance theory, we established that MIE 





is globally asymptotically stable, if Ro > 1. Numerical simulations are presented to 


confirm the analytical solutions. 


1. Introduction 


Mathematical models play a remarkable role in undestanding the dynamics of infectious 
diseases and suggest the control strategies. In the study of dynamical systems such as 
epidemiological models, the main focus is not on finding detailed solutions, but to inves- 
tigate some characteristics of the system such as existence and stability of equilibrium 


points (Lungu et al., 2007). A vector x* is an equilibrium point of a dynamical system 


x = f(x,t) if f(x*,t)=0, Vt» 0. 


An equilibrium x* is said to be stable if an arbitrary point xg of the system that starts 
near x = x* remains near it, and unstable if xo moves away from x*. An equilibrium 
is said to be locally stable if for all initial values, xo that are in a neighborhood A (x*) 
of x*, solution of the system remain near x* for all values of t. x* is said to be globally 


stable, if it is stable for all initial values xg € R”. 


Moreover, x* is asymptotically stable if it is stable and for an arbitrary initial value 
Xo, the solution of the system converges to x* as time tends to infinity. It is locally 
asymptotically stable if it is locally stable and all solutions that start in neighborhood 
of x* converge to x* as t — oo. The x* is globally asymptotically stable, if it is globally 
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stable and for all initial values xo € IR", the solution of the system tends to x* as t — oo. 
Investigation of local stability is simpler than that of global stability, though the latter is 
more preferable (Cull, 1981). Stability of system is basic requirement for its applicability 
in real-life settings, since stability justify the convergence of solutions of system towards 
a particular equilibrium point of the system (Chen, 2004). This tells us how the system 


behaves if it a solution started relatively near, but not exactly at equilibrium point. 


A number of techniques have been proposed in investigation of stability of equilibrium 
points of epidemiological models (Mpeshe et al., 2014). Linearization (Mpeshe et al., 2014; 
Li et al., 2011; Tumwiine et al., 2007a) is used on proving local stability, and Metzler 





matrix theory is used for global stability of disease free equilibrium (Mpeshe et al., 2014; 
Wang and Liao, 2012; Dumont et al., 2008; Kamgang and Sallet, 2008). Lyapunov fuctions 
has been useful tool on the study of global stability of endemic equilibrium (Kajiwara 
et al., 2015; Korobeinikov and Maini, 2004). Morever, some models are complex in such 
a way that existence and stability of equilibria cannot be investigated explicitly. Instead 


numerical simulations have been used to facilitate the purpose (Chiyaka et al., 2008). 


In this study, we investigated the existance and stability equilibrium points of mathe- 
matical model for the in-human host and in-mosquito dynamics of malaria parasites. We 
applied linearization technique to establish the local stability of MFE. We used Metzler 
theory to establish global stability of MFE. Global stability of MIE is established us- 
ing Lyapunov function in combination with LaSalle’s invariance principle. Moreover, we 


performed numerical simulations to prove the existence and stability of MIE. 


lhis paper is organized as follows: Mathematical description and formulation of the 
model for in-human host and in-mosquito dynamics of malaria parasites is presented in 
Section 2. In Section 3, we present analysis of the formulated model, whereby existence 
and stability of equilibrium points are discussed. In Section 4, numerical simulations are 
presented to prove the analytical solutions. Lastly, conclusion and direction for the future 


works is presented in Section 5. 
2. Model Formulation 


The model formulated in this study describes the in-human host and in-mosquito dy- 
namics of malaria parasites. During the blood meal infected mosquito injects sporozoites 
Sn, into human at rate abv, which then attack hepatocytic liver cells (HLCs) H, at a 
rate 615, H, and die at rate usn where £4 is infection rate of sporozoites on HLCs. The 
infected HLCs Ip, progress to hepatic-schizont T}, at a rate or, which eventually burst 








at rate 0175, to release merozoites, M. The released merozoites, then attack the healthy 
red blood cells (RBCs), R at rate 55 RM, where 05 is infection rate of merozoites on 
RBCs. The parasitized RBCs I,, progress to erythrocytic-schizonts T;., at rate oT... The 
T, burst to release either new merozoites at rate pdoT} that attack other healthy RBCs or 


gametocytes Gp, at a rate (1—p)d97;,. When uninfected mosquito bites an infected human 
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ingests the gametocytes, Gy, which develop into gametes, Gm, at the rate pquwGy, where 
p is number of bites a mosquito made during its lifetime, w is number of gametocytes in- 
gested per bite and q is probability that a mosquito bite is infective to mosquito while G, 
is number of gametocytes in blood stream. In the mosquito’s midgut the microgametes 
fuse with macrogametes, to develop into Oocysts C. at a rate a3G,,. Then, C burst to 
release sporozoites Sm, at a rate ôC, which migrates to salivary glands ready to infect a 
new host. Death rates for H, I, and Tj, are un, Lin and py, respectively. Death rates of 
R, L, and T, are Ur, Uir and Hir respectively. Sh, Sm and M die at rates Ush, Hsm and Um 
respectively. The HLCs and RBCs are recruited from bone marrow at constant rates Ap 


and A, respectively. The variables of the model are presented in Table 1. 


' Table 1. List of state variables 


Variable Descrption 


Sr: number of sporozoites in human 
H : number of uninfected HLCs 
I number of infected HLCs 
Ti number of liver schizonts 

T, number of blood schizonts 
M : number of merozoites 

R number of uninfected RBCs 
I, number of infected RBCs 
Gp: number of gametocytes 

Gm š number of gametes 

C number of Oocysts 


Un 
D yi 


number of sporozoites in mosquito 


In formulation of this model, we make the following assumptions. A cycle starts by 
a bite of infected mosquito onto uninfected human and we neglect a bite of infected 
mosquito onto an infected human host. The HLCs and RBOs recruited at constant rates 
from bone marrow and they are infected depending on their densities. Mosquito-human 
infection is does not depend sporozoites' density in salivary gland, while human-mosquito 
infection dependent of density of gametocytes in blood stream (Da et al., 2015). We also 
assume that death rates of infected cells is higher than that of uninfected ones. Also, it 


has been assumed that each of injected sporozoite and released merozoite either die or 





successfully infect HLCs and RBCs respectively. Similarly, ingested gametocytes either 


die or successfully fuse. 


Within each replication in erythrocytic cycle a constant proportion of asexual parasites 
switches to gametocytes. Finally, we assumed that existence of mosquito depends on 
human blood to develop their eggs. Based on the dynamics and assumptions stated 
above, the in-human host and in-mosquito dynamics on malaria are presented in Figure 
1. 
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Figure 1. Model compartmental diagram for in-human host and in-mosquito 


dynamics of malaria parasites 
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The detailed biological descriptions of parameters are as presented in Table 2. 


Table 2. Parameters estimates for the model (1a)-(11) 








Parameter Description Value Reference 
a: probability that a bite infects human 0.75 (Tumwiine et al., 20070) 
b: number of mosquito bites per individual 15day ! Estimated 
v : number of sporozoites injected per bite 10 — 20 (Nelson and Williams, 2014) 
£1: infection rate of HLCs by sporozoites 0.001 ulcell !day ! Estimated 
ri: number of merozoites per liver schizont 10000 (Tumwiine et al., 2014) 
Qı : progression rate of infected HCLs to schizonts 0.125 day! Estimated 
0, : rupture rate of liver schizonts 0.0975 day! Estimated 
An : the recruitmet rate of HLCs 3000 cellsday ul 1 Estimated 
un : natural death rate of uninfected HLCs 0.94 day Estimated 
Ln : death rate of infected HLCs 0.95 day! Estimated 
Ltn : death rate of liver-schizonts 0.029 day! Estimated 
G2: infection rate of RBCs by merozoites 2x10 9 ulcell !day ! ^ Estimated 
ôs : rupture rate of blood schizonts 0.115 day! Estimated 
Q : progression rate of infected RBCs to schizonts 0.145 day! Estimated 
ro: number of merozoites per blood schizont 16 (Dube et al., 2010) 
q: probability that a bite is infectious to mosquito 0.09 (Agusto et al., 2012) 
w: number of gametocytes ingested per bite 10 Estimated 
p: number of bites made by mosquito in its lifetime 3 Estimated 
A, : the recruitmet rate of RBCs 4.15 x 10! cellsul !day ! (Li et al., 2011) 
jt, : natural death rate of uninfected RBCs 0.02 day! (Dube et al., 2010) 
Hir : total death rate of uninfected RBCs 0.025 day! (Diebner et al., 2000) 
Utr : death rate of blood-schizonts 0.185 Estimated 
[lm : death rate of merozoites 48 day! (Li et al., 2011) 
Hg» : death rate of gametocytes in bloodstream 6.25 x 10 ?day ! Estimated 
63: rupture rate of Oocysts 0.05 day! Estimated 
ra : number of sporozoites per Oocyst 1000 (Nelson and Williams, 2014) 
O3: progresion rate of gametes to Oocysts 0.07 day! Estimated 
Hm : death rate of gametes in mosquito's midgut 0.052 day! Estimated 
[lc : death rate of Oocysts 0.024 day! Estimated 
Usm : death rate of sporozoites in mosqouito 40 day ! Estimated 
Hsn : death rate of sporozoites in human liver 12x10-" day"! Estimated 


p: proportion of asexual that differentiate to merozoites 0.926 Estimated 
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From the compartmental diagram in Figure 1, the dynamics of the entire in-human and 


in-mosquito malaria cycle we derive the model which is governed by the following set of 





nonlinear ordinary differential equations: 


dH 








(1a) s =A; — iS, H — ui H 

(1b) cs mh SH — (m + ta) 

(10 TE onla = (51 + na)T, 

(1d) 2 =7 017}, + pr3ó31;,. — Bo RM — ju M 
(1e) c LA, — BRM — pR 

(if) ^ —&RM — (o + bir) I, 

(1g) < =l, — (Ë + pe )T, 

(1h) T —(1 — p)rəóyT, — (qu + ug) Ch 
(1i) Cn =pquGy — (o + Hgm)Gm 

aj Co. =esG,, — (8s + He) 

(1k) sei miss = (Qp irs 

a) TR aby — &S.H — aS, 


3. Analysis of the Model 


In this section, we study the basic properties of the model system (1a)-(11). For epidemi- 
ological implications, we prove that the model system (1a)-(11) has the solution that is 


mathematically and biologically well-posed in the feasible region 


A 
Q= [a In, Th, M, R, In, Tr, Go, Gm, C, Sm; Sh) € RY : Np (t) < max wO) =) 
Hı 
A 1 A A 
N,(t) < max [N0 =} , M(t) < max [uro — k + prada | | 
H2 m H1 H2 


Gy (t) < max fao (1 — prot d , Nu (t) < max Iva, ES É =p 








Ó Ñ; 
Sm < max 500). £ n a - prada | , Sa (t) < max [s0 241 
DRIN H2 Hsh 
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where 


Hı miniga, fin +61}, fo = mint ur, Hir + do}, pa = MIN{ lgm, Hc], 


Ny(t) = H(t) + In(t) + Tht), N-t) = R(t) + L(t) + T(t), Nm(t) = Gm(t) + C (t), 


and R^ is non-negative orthant of R'*. 
Conditions for existence and stability of equilibra of the model in this region are discussed 


in next subsection. 
3.1. Existence and Stability of Equilibra 


In absence of infection, we obtain one equilibrium termed as malaria free equilibrium 


(MFE), 


A A, 
F? = (75.0,0,0, 7, 0,0,0,0,0,0,0,0) 
Hh Hr 
The stability of MFE is discussed in next subsection. 
3.1.1. Local and Global Stability of MFE 


We establish the local stability of E? by investigating the signs of the real parts of the 
eigenvalues of the Jacobian matrix of the system at E?. Jacobian matrix of system (1a)- 
(11) at E? is given by 


— Ü 0 0 0 0 0 0 0 0 0 — zı 
0 -z2 0 0 0 0 0 0 0 0 0 AI 
0 O1 —23 0 0 0 0 0 0 0 0 0 
0 0 710, —z4 0 0 prodg 0 0 0 0 0 
0 0 0 —z; —ur 0 0 0 0 0 0 0 
J(E®) = 0 0 0 Z5 0 —z 0 0 0 0 0 0 
0 0 0 0 0 O9 —27 0 0 0 0 0 
0 0 0 0 0 0 28 —z 0 0 0 0 
Ü) 0 Ü) Ü Ü) 0 0  pqw —zo 0 Ü) Ü) 
0 0 0 0 0 0 0 0 O3 —2, Ü 0 
0 0 0 0 0 0 0 0 0 T303 22 Ü 
0 0 0 0 0 0 0 0 0 0 0 —243 
where 


Ap Ns Ny 
A= — | ES as SL es. 2S bm. SS 

Hh Hr Hr 
Zg =A + a 27 = 02+ fer, %= (1 E p)r202, Zo = QW + Hgb, 210 = Q3 + Hgm, 


A 
(2) 24, =63+ Me, 212 = GU + Hsm, 213 = A + Hsh 
h 
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The MFE is locally assympotically stable if and only if trace of J( E") is strictly negative 
and determinat of J(Eo) is strictly positive. We obtain the following results, 


trace(J(E°)) = — [(un + Hr) + (uan + Htr) + (Bh + Hir) + (a1 + 02 + 03) + (01 + Ó> + 63) 


Bi An a BoA, 


(3) + (Um + Hgb + Be + Hsm + Tuas) + 
Hh Hr 








+ (qu + pgm + +av)] < 0 


and 


det(J(E°)) =(GiAn + pan pin) (01 + pan) (Gi + pin) (qoo + Hgo) (ors + gm) (0s + pic) (av + Hsm) 
(4) (Hir + 02) (02 + pir) (P9 + Hmhr) |l — Ro] > 0 
where 


BoA, O2 pr202 


Riss > 
i BoA, a HmhHr (a T liir) (55 T Ltr) 


Equation (4) holds only if Ro < 1 ; and because of this requirement Ro is interpreted as 
the basic reproduction number. This leads us to the following theorem. 
Theorem 1. The malaria-free equilibrium, E? is locally asymptotically stable when Ro < 


1 and unstable otherwise. 


3.1.2. Global Stability of MFE 





We establish the global stability of E? using the Metzler matrix theory technique used in 
Castillo- Chávez et al. (2002); Kamgang and Sallet (2008); Mpeshe et al. (2014). In this 
approach, we re-write the model system in the form: 


afm = Ai(z)(X, — Xeon) + Arlt) Xe 





A = A» (x) X, 





where X, is the vector of uninfected classes and X, is the vector of infected classes. For 


our case, we have 


(5) Ar LI (H, R) and Xe — Uis Thy M, I. Gp, C gi C, Du Sh) 
Ny Ny 

6) Xma = (2, =] 
Hh Hr 

and 
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0 0 0 000000 -6H 
(8) Anla) = * 
00 -6R 00 0 0 0 0 0 
and 
—w, 0 0 0 0 0 0 0 0 6H 
Qy —W 0 0 0 0 0 0 0 0 
0 T104 — W83 0 pr30» 0 0 0 0 0 
0 0 BR —w, 0 0 0 0 0 0 
0 0 0 — 0 0 0 0 0 
(9) Ad(x) = “< 
0 0 0 0 We —u; 0 0 0 0 
Ü) Ü) 0 Ü) 0 pqw —wg 0 Ü) 0 
0 0 0 0 0 0 O3 —W9 0 0 
0 0 0 0 0 0 0  Taós —149 0 
0 0 0 0 0 0 0 0 W11 —Wi1e 
where 
En _ _ BoA, o 
Wy —0O + Hih, W2 = 01 + Mtn, W3 = + Um, W4 = O2» F lirs 
Us =01 + Mir, We = (1 — p)r2ô2, W7 = qu + Hgb, Wg = Q3 + Hgm, 
abv 
(10) Wg =03 + Ue, Wino = QV + Hsm, W11 = —— and wiz = Pí H + pan 


Sm 


It can easily be seen from (7) that, all eigenvalues of A; are real and negative. 


So, the system 





= Ai(x)(X4 — X rin) + Ai(x) X. 


is globally assymptotically stable at Xgo. From (9) and (10) it can be observed that all 
off diagonal elements of Ay are non-negative. Therefore, Ay is a Metzler stable matrix. 
Thus, the MFE is GAS. To investigate under which conditions MFE is GAS, we need to 


prove the following proposition . 


Proposition 1. (Kamgang and Sallet, 2008; Dumont et al., 2008) 

A B 
Let M be a square block decomposed Metzler matriz: M = OD with A and D 
square matrices. Then M is Metzler stable if and only if matrices A and D — C ATHA are 


Metzler stable. 
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For this case we have M = Ag, and A = 0 yoy —us Ü  pra3ó9 


0000 &H 0 0 0 0 "wg 
0000 0 0000 0 
p=] 0000 O0 ; C=] 0 000 0 
0000 0 0000 0 
0000 0 0000 0 
—u,; 0 Ü) 0 0 
pqu —wg 0 0 0 
and D = 0 Q3 —wWg9 0 0 
0 Ü | raOQ3 -—1190 0 
0 0 0 Wir | —12 
immediately we have 
E 0 0 0 10101026, Ar B4 Bo 
W1Woup(W3W4W5 Ur — D», pra0203) 
"o Wi? . —1A43 0 0 0 
i 0 Q3 —Wi14 0 0 
0 0 Win —Wi¢ 0 
0 0 0 W17 — W118 
Definition: 


A Metzler matrix M is said to be stable if all of its diagonal elements are negative. 


By that definition, A is Metzler stable matrix, and D — CA !B is Metzler stable matrix 
if and only if 





Ó Ap, 
(11) r1010302w6A 5 A, p1 2 +0 
W1 Woh (W3W4W5 Ur — Da, praó203) 
which holds only when 
A, pr3ó 
(12) B» DT 2020/2 PM 
W3W4W5 fUr 
Using ws, w4, and ws as given in equation (10), we get 
B2 Nr 
Oi3pr20» 
(13) pr e 


(= + T (aa Mir) 0t Hi) 
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equivalently, 
Ro < 1. 


This leads us to the tollowing theorem. 


Theorem 2. The MFE of the model system (1a)-(11) ¿s globally assymptotically stable 
in Ql if Ro < 1 and unstable if Ro > 1. 


3.2. Existence of MIE 


The model has one positive malaria infection equilibrium E* which is given by 


E = AH ob MV Oso uo S.) 


where 
"d Ap "- By ApS), El B1 Apna, SF, 
H* = BN eae an ERG In ZU PO TOU EU e PRONTO SCENES. 
Diog + Hn (ar + Min) (P155 + un) (01 + pan) (0a + Bin) O15, + un) 
: Ny š BoNoM* a Tu (Bəs A|, + Hrm) M* 


7 Boa M* + u ” i (ao + pir) (Go2M* + ur) ” praó»(B5 M* + ur) 
x __ (1 = p)Ro (BoA, 25 Hy] ) M* I pqu (1 7 p)Ro (BoA, a lir Hm) M * 
Gt = ——— —— Ga = | —_ | | —— |, 
p(qu- Hg) | (82M* + ur) Q3 + Hgm] |p(qu + Hg)  (B2M* + ur) 
PR: | pq | | (1 — p)Ro i ei] 
63 + ue [3+ Hgm] LP(qo + ugo) (BaM* +u) J 


o "sôs as | pqw | | (1 — p)Ro cem 
™ QV + Hsm Ó3 + Hc O3 + Hgm] |p(qw t+ ug) (6oM* + n.) E 


S% and M* are positive solutions of F(5S7) = 0 and G(M*) = 0 respectively, where 
F(S:) = AsS2 + AgSt + A1, G(M*) = B3M* + B4M* + B, and 


A3 = Piptsn, Ag = B1(An — aby) + usns, A1 = —abvypn, 
Bs = Ban, Ba = ~[(Ba Ida = 1) + 2X], Bi = =p" and 
\(S*) = — 
Si) Bi Sp + un 
Now we need to determine necessary and sufficient conditions for existence of malaria- 


infection equilibrium E* by proving the following theorem 

Theorem 3. The model system (1a)-1l has a unique malaria infection equilibrium 
EH P cuu T Gu 15 ) 

if Rog» 1, An > abv, Aa3S;7 + A25% + A, = 0 and B4M** + B3M* + B4 = 0 have roots 

55 > 0 and M* > 0 respectively. 
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Proof: 
Let E* = (H* I, TE, M*, R*, Pr T, Gu, G>, C*, S*,, Sr) be malaria infection equilibrium 
of the system (1a)-(11). Substituting the expression for H* into equation (11), we have 

Bios n 


aby — —— = 
Di Sr + Hh 


— LshOn, =) 
which yields to 


(14) Din S, + [P (A, — abv) + nsn] Sp — abv, = 0 


Since A; < 0 and A; > 0, then the quadratic equation (14) has unique positive root S% 


given by 
(15) Ta _ [83 (A m abv) zi Hsnn| 25 VA, 
ý 2 Bihan 
where 
(16) A = (81(An — abv) + Herth) + 481abvusn n 


only if A» > 0. 


Hence, Ag = ñi (A, — abv) + uan, > 0 only if A, > aby (recruitment rate of uninfected 


HLCs is greater than recruitment of sporozoites into human liver). 


Substituting expressions for T7, TZ and R* into equation (1d) gives 


(17) 
Di As 0371015; Ro(Ba X, + Umr) M* B3. M* " 
(01 + Lun) (o3 + un + dn) (Gi Sh + pn) (B5.M* + ur) BoM* + br 
; P1 Ap 04T101 I 
Letting Ü = —— — PP equation (17) becomes 
: (01 + un) (ai + pon + dp) 9 QT) 
(18) 055 [Rol BoA, T Hals) v Boh, = Hml D> M * RUE Hr )] M* —() 
(155 + un) (B5 M* + ur) 
which can be further simplified to 
(19) Bau, M (82. sp Inr) (Ro 1) E Bod") M se [pA = 0 
where : 
(P157 + un) 


Equation (19) has a unique positive real root M* given by 


x (BoA, E Hoe US 1) E B3A*| RE V/A» 
m E 2 Bo bm us A* 
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where 
Ao m (Ro S 1)" (B3 A, SP TNR T 252A" | BoA, RE rape Ue _ 1)| "E BA Es 4 bo Lm [br AT 


only if Ro > 1. 


This is condition for the existence of malaria infection equilibrium. 
Therefore, if Ro > 1, A, > aby and quadratic equations A357? + A557 + A; = 0 and 
B3M* + B3M* + B, = 0 have respectively positive roots S* and M*, with 


ot — v fen (Ap, ES abv) T js dus] T VA, and M* = (89, 3b Lmbr)(Ro m 1) "e B3A*| NE V A2 
I 201 [sh 2 D2 [bm [br A 


3.3. Global Stability of MIE 


When Ro < 1, MFE is locally asymptotically stable. This suggests local stability of the 
MIE when Z > 1 (Van den Driessche and Watmough, 2002). 


50 we only investigate the global stability of the MIE. We adopted the techiniques used 
by Pedro et al. (2014). 


Theorem 4. if Ro > 1, the model described by equations (1a)-(1l) has unique positive 
MIE, E*, such that 




















aes Jor O < S, < S, and 0 < H < H* 
EE ET for0< S, < S, and 0 < H < H° 
zu il for0<M<M* and 0 < R < PR 
ree el for0<M<M* and 0 < R < R° 


Then, E* is globally asymptotic stable in Qe. 





Proof: To estabilish the global stability of MIE, E*, we define the Lyapunov function of 


the form 


Le) = X z (zi — 2; In j lor vex. Deel 


as proposed by Castillo-Chávez et al. (2002) where z; is a number of cells in the i‘ class, 
x; are equilibrium values and z; are constants. This approach has been found useful for 


more complex compartmental models of in vivo dynamics (Korobeinikov, 2004). 
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Now, we constructed the following lyapunov function 
H I T, 
L=xy (u- H* ln Z) + 2 (i - I;ln E) + 23 (n - T; In rJ 
M R I, 
=) + 25 (R-P) + 26 Ú - mg) 


l; * Gi * Gm 
Tor (r. — T* In 2 + 2g (a, - G; In =) + £9 (a. — G; ln 2) 





+24 (m-ar 














C Pu S 
(21) +210 (c - C" In =) + Zil (s, m D In Sm) + 219 (s m s In =) 
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By hypothesis of Theorem 4, we have 




















a ot H* M* | | R*M* 
R* | | R*M* S; | | S;H* 


where equality applies only when H = H*, S, = Sr, M = M*, R = R* Therefore 
f(Q) < 0 for all H = H*, Sn = Sr, M = M*, R = R*. 
dL dL 
Hence, — < 0 for all (H, I4, T4, M, R, Ip, Te, Gb, Gm, C, Sm, M > 0 and — = 0 only 
when T = I” Pl phi Ph Pb, M = M' h= Krd = Tr Guo Gmn = Gr 
C = C*, Sm = 55, Sh = Sy, and E* is the only equilibrium TN "i the system on this 
d.L 
plane. Therefore, the largest compact invariant set in Q such that — = 0 is the singleton 


LE* Y which is the MIE. LaSalles invariant principle (LaSalle, 1976) guarantees that E* 
is globally asymptotically stable (GAS) in the interior Q of Q. 


4. Numerical simulations 
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In this section, we perform some numerical simulations of the model (1a)-(11), to illustrate 
the dynamics of model using MATLAB symbolic package run in intel (R) Pentium (R) 
CPU B980 2.40GHz, 2.40GHz, 4.00GB machine. The initial values used in simulation 
of this model are largely assummed to allow computer executions, and their values are 
listed in Table 3. 


Table 3. Initial values of variables of the model (1a)-(11) 


Variable H I, | T, | M R 1, T. G Gm |C Om, | S 
Initial values | 3000 0 |0 | 2000 | 500000 | 0 | 1000 | 3000 | 1500 | 1000 | 2000 | 2000 

















Although the decision on values of parameters for the in vivo dynamincs is challenging 
(Chiyaka et al., 2008), the numerical values of parameters used in the numerical simula- 
tion of this model are presented in Table 2. These values are either estimated or taken 


from various articles among existing literature. The reason why some parameters values 





have been estimated is that modelling of liver and mosquito stages dynamics of malaria 


parasite have not been done or the parameter values found in existing literature are not 





suitable in our model. Even those that have been taken from other related studies may 
not be as accurate as we need for our mathematical forecasts. However, the main issue 
here is the effect of these parameters on the basic reproduction number, which gives the 


clues on how to eradicate or control the disease (Chiyaka et al., 2008). 
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Figure 2. Variation of populations at Exo-erythrocytic cycle with time 


It is observed from Figure 2 that population density for sporozoites, uninfected HLCs, 


and infected HLCs vary with time and attain constant values (malaria infection point). 





However, sporozoites injected into human starts by falling within very short time before 
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starting to rise. ‘This fall is probably due the fact that when sporozoites injected into the 
human, they travel to the liver through bloodstream where they ingested by phagocytes 
(Smyth and Wakelin, 1994) or they probably die due to change of environment from 
mosquito’s salivary gland to human bloodstream. ‘Then its population increases after they 
succefully reach the liver and start the asexual replication (exo-erythrocytic schizogony) 
within HLCs. In contrary to population of sporozoites, population of uninfected HLCs 
decreases with time until it reaches its equilibrium value. ‘This population decreases 
because of infection of HLCs by sporozoites, which on other hand cause rise in population 


of infected HLCs and liver schizonts. 
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Figure 3. Variation of populations at Erythrocytic and Sporogonic cycles with time 





Figure 3b shows that the density of Gametes rises to its maximum value before falling 
to its equilibrium value, while those of oocysts and in-mosquito sporozoites increase with 
time until they attain constant values (values at malaria infection equilibrium). The 
decrease of gametes may be it is due to the formation of ookinetes. Therefore from the 
Figure 3, we conclude that the malaria-infection equilibrium, E* for this model exists. 


Now let us assess for stability of E*. 


Using the parameters values given in Table 2, we obtained Rp = 1.59025 > 1. Thus, by 
Theorem 4 implies that the malaria infection equilibrium £* is globally assymptotically 
stable as depicted in Figures 4, 5 and 6. It has observed that with different initial values, 


solutions trajectories for all state variables converge to malaria infection equilibrium. 
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Figure 4. Numerical simulation to show global stability of MIE for variables in 


exo-erythrocytic cycle 
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Figure 5. Numerical simulation to show global stability of MIE for variables in 


erythrocytic cycle 
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Figure 6. Numerical simulation to show global stability of MIE for variables in 


Sporogonic cycle 


5. Discussion and Conclusion 


The mathematical model for in-human host and in-mosquito dynamics of malaria para- 
sites was developed and analyzed. The model involved three main phases in life cycle of 


malaria parasites. We considered four, five and three compartments in the liver, blood 





and mosquito stages respectively. 


In analysis of the model, two steady states, malaria-free equilibrium (MFE) and malaria- 





infection equilibrium (MIE) were determined. The threshold, Ro, was obtained and found 
to a function that depends only on parameters in erythrocytic phase. This implies that 
the erythrocytic invasion may propagate without new infection from the liver (implying 
that even when an individual is not bitten by the mosquitoes, s/he may maintain some 


level of malaria in the blood). 


20 MOHAMED A. SELEMANI, LIVINGSTONE $. LUBOOBI AND YAW NKANSAH-GYEKYE 


Stability (in terms of Ro) of equilibrium points was established. The necessary condition 
for stability of MFE was established using trace-determinant of jacobian matrix of the 
model evaluated at this point showed that, MFE is locally assymptotically stable provided 
Ro < 1 and unstable otherwise. The global stability of this equilibrium was investigated 
using Metlzer matrix technique, and proven that MFE is globally asymptotically stable 
when Rp. < 1. MIE exists only if the recruitnent rate of sporozoites into human host 
less than recruitment hepatocytes liver cells and Ro > 1. Global stability of this was 


investigated using Lyaponuv function. 


An insight of dynamics of malaria parasites within human host and within mosquito is 
significant in development and assessment of transmission blocking intervetions (T BIs). 
Merozoites play an important role in propagation of malaria infection in human, and 
they initially procuded in the HLCs after invasion of sporozoites. This may suggest 
that blocking this invasion to be one of the best targets for TBIs as it will significantly 
inhibits the infection of HLCs, and eventually the production of merozoites from the liver 
schizonts. Therefore, it reduces the possibility for infection of RBCs by merozites from 
the liver. However, as it has been stated earlier that infection of RBCs by merozoites may 
propagate without a new infection from the liver, but this would occur only when initial 
invasion of RBCs by merozoites from the liver was successful. Therefore, implementing 
the TBIs at liver stage will probably reduce possibility of having erythrocyte invasion of 


merozoites and finally the human-mosquito transmission may be stopped. 


This work provides a basic model for studying the in-human host and in-mosquito vector 
dynamics of malaria parasite. At this time where malaria eradication is on world agenda, 
this work may be used as starting point to examine how and which are new control 
strategies of malaria can be established to overcome the disease. It will be useful to 


study the effect of immune response and/or treatment in the extension of this model. 
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